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ABSTRACT 

We present an overview of and first results from the OMEGA survey: the OSIRIS 
Mapping of Emission-line Galaxies in the multi-cluster system A901/2. The ultimate 
goal of this project is to study star formation and AGN activity across a broad range 
of environments at a single redshift. Using the tuneable-filter mode of the OSIRIS 
instrument on GTC, we target Ra and [Nil] emission lines over a ^0.5 x 0.5 deg^ 
region containing the 2 : ^ 0.167 multi-cluster system A901/2. In this paper we describe 
the design of the survey, the observations and the data analysis techniques developed. 
We then present early results from two OSIRIS pointings centred on the cores of the 
A901a and A902 clusters. AGN and star-forming (SE) objects are identified using the 
[Nil] / Ra vs. Wjia (WHAN) diagnostic diagram. The AGN hosts are brighter, more 
massive, and possess earlier-type morphologies than SE galaxies. Both populations 
tend to be located towards the outskirts of the high density regions we study. The 
typical Ra luminosity of these sources is significantly lower than that of field galaxies 
at similar redshifts, but greater than that found for A1689, a rich cluster at 2 ; ^ 0.2. 
The Ra luminosities of our objects translate into star-formation rates (SERs) between 
^ 0.02 and 6 Mq yr“^. Comparing the relationship between stellar mass and Ha- 
derived SER with that found in the field indicates a suppression of star formation 
in the cores of the clusters. These findings agree with previous investigations of this 
multi-cluster structure, based on other star formation indicators, and demonstrate the 
power of tuneable filters for this kind of study. 

Key words: Galaxies: clusters — Galaxies: AGN — Galaxies: star formation. 


1 INTRODUCTION 

Galaxies are not randomly distributed in the Universe: while 
most are part of groups and clusters, they are also found in 
isolation or in the field. The observed properties of galaxies 


* The first two authors have contributed equally to this paper, 
t E-mail: ana.chies.santos@iag.usp.br 


are strongly dependent on the environment in which they 
reside ( Dressler|1980 ). In the local universe, galaxies in low- 
density environments are more likely to have spiral mor¬ 
phologies, blue colours and be forming stars, whereas those 
living in high-density regions are largely bulge-dominated, 
red and passive ( Balogh et al.|[2004 Bamford et al.|[2009 ). 
These two populations may be connected by one or more 
evolutionary processes by which galaxies falling into denser 
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environments have their star formation quenched and their 
internal velocity dispersions increased. 

Both internal and external processes can drive this 
transformation. External processes such as ram-pressure 
stripping (Gunn & Gott 1972 Abadi et al.||1999 Bekki et 


al.|2002 ), tidal interactions ( Larson et al.|1980| ' Bekki[1999 )^ 


galaxy harassment ( Moore et al.||1996 ), or mergers (Barnes 
1992 Bekki||1999 ) can produce such morphological changes 


and cause the removal and/or consumption of gas, leading 
to the halt of star formation. On the other hand, processes 


related to internal galaxy properties, such as mass ( Haines et 
al.|2006 ) and mass-dependent feedback mechanisms (Bundy 


et al. 2006) can also be responsible for the quenching of 


star formation, as can feedback from super novae and ac¬ 
tive galactic nuclei (AGN; Booth Schaye|2009 Newton & 


Kay|2013 ). In fact, environment also seems to be connected 


with the probability of a galaxy to host an AGN: while low- 
luminosity AGN do not show particular preference for any 
kind of environment ( Martini et al.|2002 Miller et al.|2003' 


Martini et al.|2006 ), luminous AGN are preferentially found 


in the field (Kauffmann et al. 2004) and groups (Popesso 
Bavian^|2006 ) rather than in clusters. Peng et al. (2010) 


and Thomas et al. (2010) showed that both stellar-mass and 


environmentally driven quenching of star formation are im¬ 
portant in governing the evolution of a galaxy and that their 
effects can be separated from each other. 


Vogt et al. (2004) and Jaffe et al. (2011) have shown that 


galaxies falling into clusters experience a removal of their 
gas, which in principle should lead to a gradual decline of the 
star formation activity. However, the average star formation 
rate (SFR) in star-forming (SF) galaxies remains roughly 
constant with environment, and it is only the fraction of star¬ 
forming ga laxies that changes ([Balogh et aL||2004 Verdugo 


et al.|200^ Poggianti et al.|2008 Bamford et al.|2008 ). The 

extent of ongoing star formation within a galaxy has also 


been shown to correlate with environment (Koopmann V 


Kenney||20Q4| Moss V Whittle||2000| Bamford et al.||2007 ), 
implying a truncation of the star formation taking place first 
in the outer parts of the infalling galaxies. In fact, the last 
episode of star formation in cluster disc galaxies is seen to 


be concentrated towards the central regions ([Rodriguez Del 
Pino et al. ]|2014| [Johnston et al.||2014| ). 


In order to understand the full degree of transformation 
and the relative role played by external and internal pro¬ 
cesses, we selected the field of the A901/2 multi-cluster sys¬ 
tem at z ~ 0.167. This large structure has been the subject 
of study of the STAGES project ( [Gray et al.|2009[ ), where a 
wealth of data exists: an 80-orbit F606W HST/ACS mosaic 
covering the full ~0.5 x 0.5 deg^ (~5 x 5 Mpc^) span of the 
system, complemented by extensive multi-wave length obser¬ 
vations with XMM-Newton, GALEX, Spitzer, 2dF, GMRT, 
Magellan, and the 17-band GOMBO-17 photometric redshift 
survey ( Wolf et al.|2003 Wolf et al.|2004 ). The broad range 
of galaxy environments and luminosities sampled by these 
observations provide the perfect laboratory to study both 
obscured and unobscured star formation, stellar masses, 
AGN activity, and galaxy morphologies. 

Previous work in the A901/2 field has included studies 
of both star formation ([Gray et al.[[2004) and AGN activ¬ 
ity and how they relate to galaxy environment. [Wolf et al.[ 
(2005) discovered red spirals in the system which suggested 


galaxies. Gallazzi et al. ( [2009 ) found an overall suppression 
in the fraction of star-forming galaxies with density by com¬ 
bining the UV/optical SEDs from GOMBO-17 and Spitzer 
24 pm photometry. The surviving star formation was shown 
to be obscured and hosted by a relatively large fraction of red 
spirals located primarily in the cluster infall regions ( [Wolf 
et al.|2009 ). In a more recent work, Bosch et al. ( 2013| find 
that this suppression of star formation is partially driven by 
ram-pressure stripping, which could lead to the production 
of red spirals and, ultimately, SOs. Although the integrated 
star-format ion properties have yielded important results, the 
actual spatial scale of the star-forming regions and how it 
is affected by the environment remains unknown. Likewise, 
initial attempts at quantifying the AGN activity in the field 


of A901/2 were made by Gilmour et al. (2007) and Gallazzi 
[et al.[ ( [2009[ ) using XMM-Newton X-ray and Spitzer 24 pm 
data, respectively. However, the lack of emission-line diag¬ 
nostics means that the AGN census remains incomplete. 

In order to obtain spatially-resolved emission-line di¬ 
agnostics, we have used the Gran Telescopic Ganarias 
(GTG) Optical System for Imaging and low-Intermediate- 
Resolution Integrated Spectroscopy (OSIRIS) in tuneable 
filter mode to obtain very deep spatially-resolved emission¬ 
line images of the complete A901/2 STAGES field. We have 
targeted the Ha (Ao = 6563 A) and [Nil] (Ao = 6583 A) 
lines. The former is the best optical tracer of SFR ( Kennicutt[ 
1998), whereas the combination of the two lines provides one 


of the most reliable optical AGN diagnostics available ( [Gid 
Fernandes et al.[[2010 ). 


reduced current star formation rates in a large fraction of the 


In this paper we present the design of the OSIRIS Map¬ 
ping of Emission-line Galaxies in A901/2 (OMEGA) survey, 
the data reduction procedures adopted, and a first analy¬ 
sis of two of the twenty pointings, demonstrating some of 
the science that the survey can achieve. OMEGA has been 
designed to address star formation and AGN activity as a 
function of galaxy properties and environment in A901/2. 
In the present work, we show a tentative census of AGN 
and star-forming galaxies in the highest-density cluster re¬ 
gions. A full environmental analysis of both integrated and 
spatially-resolved galaxy properties will be presented in sub¬ 
sequent papers. 

We adopt an Hq = 70 kms“^Mpc“^, Qm = 0.3 and 
Cl A = 0.7 cosmology throughout. 


2 DESCRIPTION OF THE PROGRAMME 

AND DESIGN OF THE OBSERVATIONS 

The OMEGA survey was designed to obtain very deep 
spatially-resolved emission-line images for the A90I/2 multi¬ 
cluster system and, in turn, recover low-resolution spectra in 
a wavelength range that covers Ha and [Nil] for the galaxies 
of the system. Moreover, OMEGA is intended to sample the 
16.0 ^ ruR ^ 23 magnitude range. This roughly corresponds 
to stellar masses of 9.0 ^ log (M*/Mq) ^ 11.5 at the clus¬ 
ter redshift, reaching lower SFRs than possible at high-z. 
The survey is based on a 90 h ESO/GTG Large Programme 
allocation (PL A. Aragon-Salamanca). 

In order to meet the objectives of OMEGA, we make use 
of the tuneable filter capabilities of the OSIRIS instrument 
( Gepa et al.|20I3a] Gepa|20l'3b ), located at the Nasmyth-B 
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focus of the 10.4 m GTC, at the Roque de los Muchachos 
Observatory on La Palma. 

Tuneable filter imaging is a method that lies between 
classical narrow-band imaging and spectroscopy, in terms 
of observing time costs and the quality of the informa¬ 
tion obtained. Although the spectral resolution is relatively 
low, the wavelength sampling is quite sufficient for reliable 
star formation and AGN activity estimates for most of the 
cluster members. Mapping the spectral region around Ha, 
[Nil] A6548 and [Nil] A6583 for all the galaxies in the A901/2 
system would be much more expensive with traditional spec¬ 
troscopy. A further advantage of tuneable filter imaging, 
with respect to slit spectroscopy, is the spatial information 
it provides. While long-slit spectroscopy normally only pro¬ 
vides information along one of the axes of the galaxies (usu¬ 
ally chosen to be the major axis), tuneable filter imaging 
provides 2D maps. The tuneable filter technique has not 
been much employed previously in this field. A heroic at¬ 
tempt, includes the GADIS survey of emission line galaxies 


measuring SFRs between z = 0.25 and z = 1.2 (Hippelein 
et al.||2003| |. 


Our increased spectral resolution over standard narrow- 
band Ha imaging surveys (see, e.g., Kennicutt|1992 Gallego 


et al. 1997| has the advantages of avoiding contamination 


from the [Nil] line in the star formation estimates and en¬ 
abling AGN to be identified ( Gid Fernandes et al.|2010 ). Fi¬ 
nally, while integral field spectrograph (IFS) surveys, such 
as GALIFA and MaNGA, also provide spatial information, 
their fields of view are much smaller than that of OSIRIS. 
Even large IFSs - such as MUSE, recently installed on the 
VLT - have a field of view that is ~50 times smaller than 
that of OSIRIS. The full unvignetted field-of-view (EOV) of 
OSIRIS is 7.8 X 7.8 arcmin^, imaged using two 2048 x 4096 
Marconi GGDs with a 9.4 arcsec gap between them. The 
plate scale, using 2x2 pixel binning, is 0.25 arcsec pixel”^. 
We can therefore obtain spectral and spatial information 
for a much larger number of galaxies, and thus gain a much 
more complete picture, than is possible with current lESs. 

The OSIRIS tuneable filter system consists of two re¬ 
flecting plates (a Eabry-Perot etalon) working in a colli¬ 
mated beam. The resulting interference allows only specific 
narrow wavelength ranges to be transmitted. One can tune 
the properties of the filter by modifying the separation be¬ 
tween the optical plates of the device. Standard narrow-band 
order-sorting filters (OSE) are used to isolate specific inter¬ 
ference orders. Reliable observations are limited to a circular 
field of view of radius 4 arcmin, within which they are as¬ 
sured to be free from contamination by other orders. OSIRIS 
has two different tuneable filters, one for the blue and one 
for the red wavelength range. We use the latter, which is 
capable of tuning filters to wavelengths in the 6510-9350 A 
range with a spectral resolution of I4A EWHM. 

A given radius from the optical centre in the image 
plane corresponds to a specific angle in the incoming col¬ 
limated beam, which is incident on the Eabry-Perot etalon. 
As the transmission wavelength of the interferometer de¬ 
pends upon incidence angle, the wavelength of the trans¬ 
mitted light varies as a function of position in the image. 
The result is a circular pattern with rings of constant effec¬ 
tive wavelength. A, which depends on the radius, r, from the 
optical centre. In order to achieve the best accuracy in the 
estimation of the wavelength we use the following expression 


( Gonzalez et al.|20I4 ): 

A = Ao — 5.04r + a3(A)r , (I) 

where 


a 3 (A) = 6.0396 - 1.5698 x I0“^A + 1.0024 x I0“^A^ (2) 


where A is in A and Ao is the effective wavelength at the 
optical centre. Since the useable EOV is limited to a radius of 
4 arcmin, the maximum variation in the effective wavelength 

is 80 A. 

In Eig. [^we show the OSIRIS observed fields as circles 
8 arcmin in diameter, overlaid on top of the footprint of the 
80-tile STAGES HST/ACS mosaic (Gray et al. 2009) and 
ground-based GOMBO-I7 R-band image ( Wolf et al.|2004 ). 
Gluster galaxies with ttir < 23.5 are also plotted. Given 
the available observing time, the initial plan was to tile the 
field-of-view uniformly in a 4 x 5 grid. However, the North- 
East and North-West tiles contain bright stars that would 
have produced reflection ghosts, seriously compromising the 
quality of the data for these tiles. Eor this reason, we decided 
to re-locate these two pointings to cover the regions with 
the highest galaxy densities, which correspond to Eields 21 
(A90Ia) and 22 (A902). These fields are referred to as E2I 
and E22 respectively. 

Fig- 0 shows the spectroscopically determined redshift 
distribution of galaxies within 1.2 Mpc diameter apertures 
centred upon the three main cluster cores. The redshifts 
were drawn from a spectroscopic survey of the 300 brightest 
cluster galaxies in the region obtained with the AAT 2dE 
spectrograph (Gray et ah, in prep.). The dashed lines indi¬ 
cate the redshift range probed by our observations. Within 
these redshift ranges we target both the Ha and [Nil] A6583 
emission lines. The z ranges are slightly narrower for the 
two pointed observations of the A90Ia and A902 regions 


(E2I and E22). 

Eor an optimal deblending of Ha and [Nil], the tuneable 
filter Eull Width at Half Maximum (EWHM) bandwidth is 
set to 14 A, and the spacing between successive wavelengths 
to 7 A. Therefore, in order to cover the desired wavelength 
range in E2I and E22 we require 14 and 12 wavelength set¬ 
tings, respectively, whereas for the rest of the fields we need 
16 wavelength settings (see Eigs. and and Table [^. 
Moreover, for each wavelength setting, 3 dithered images 
are taken typically with a ~ 40pix (10^^) separation. This 
increases our number of images (and individual wavelengths 
sampled) by a factor of three and ensures that the inter- 
GGD gap is imaged. A total of 42 and 36 images have been 
obtained for Eields 21 and 22 respectively, and 48 for the re¬ 
maining 18 fields. Note however, that the wavelength varies 
across the field, according to Eq. In Eig. we show an 
illustration of this dependence. 

Since the wavelength varies across the EOV, the central 
wavelength Ac of each wavelength setting needs to be opti¬ 
mised in such a way that the desired wavelength range is 
covered for as many galaxies as possible. The average wave¬ 
length Aave of each setting was calculated as the mean A 
within the central 5.6 arcmin diameter circle of each point¬ 
ing (the approximate area uniquely covered by one tile). Eor 
fields E2 to E20, 7609 ^ Aave (A) ^ 7714, ensuring that 
both Ha and [Nil] are observed in the 0.1594 ^ z ^ 0.I7I8 
redshift ranges. This corresponds to rest-frame velocities in 
the range —1447 ^ u(km/s) ^ 1751, where u = 0 cor- 
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Figure 1. The OSIRIS observed fields (labeled 2-3, 5-22), represented as circles of radius 4 arcmin, overlaid on top of the footprint 
of the 80-tile STAGES HST/ACS mosaic and a COMBO-17 R-band image. The very small black circles mark the positions of cluster 
galaxies with ttir < 23, where cluster membership is defined as in [Gray et al.| ( [2009| > using the COMBO-17 photometric redshifts. 


responds to z = 0.167, the average redshift of the struc¬ 
ture. For field F21, 7595 ^ Aave(A) ^ 7686, and for F22, 
7623 ^ Aave(A) ^ 7700, thus covering a slightly narrower 
redshift/velocity range for these two fields (see Fig. and 
Table for details). Note that these are average values: the 
exact wavelength range covered for each galaxy depends on 
its position with respect to the centre of the pointing. More¬ 
over, given the significant tile overlap, a sizeable fraction of 
the galaxies will be observed in more than one pointing, in¬ 
creasing the available wavelength range and the number of 
independent wavelength samples obtained. 


A summary of the observational setup is given in Ta- 
bleg The observations were executed in ‘Observing Blocks’ 
(OBs), as shown in Table and illustrated in Fig. They 
were carried out during the first halves of 2012 and 2013, 
and completed in February 2014. 


3 DATA REDUCTION 

In this section we describe the data reduction methods we 
have developed for the OMEGA data. 


3.1 Imaging Data Reduction 


The data reduction was carried out with the Osiris Ofhine 
Pipeline software (OOPs; |Ederochte||2012 ). The OOPs re¬ 
duction consists of two steps. In the first step we prepare 
the calibration files ‘masterbias’ and ‘masterflat’ combining 
the individual bias and flatfield frames taken at the tele¬ 
scope. Flatfields were obtained using dome flats with the 
Alter tuned to the same wavelength as the science obser¬ 
vations. Usually 3-5 flat frames were taken each night at 
each wavelength. In the second step we reduce the science 
frames: we remove the overscan and subtract the bias from 
the science frames and divide these frames by the master- 
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Table 1. Summary of the Observations. 


OBs (Main Survey) 

OBs (A901a/F21) 

OBs (A902/F22) 

Ac (A) 

Aave (A) 

Order Sorting Filter 

Exp. Time (s) 


1 


7615 

7595 

f754/50 

3x200 


1 


7622 

7602 

f754/50 

3x200 

1 

1 


7629 

7609 

f754/50 

3x200 

1 

1 


7636 

7616 

f754/50 

3x200 

1 

2 

1 

7643 

7623 

f754/50 

3x200 

1 

2 

1 

7650 

7630 

f754/50 

3x200 

2 

2 

1 

7657 

7637 

f754/50 

3x200 

2 

3 

1 

7664 

7644 

f754/50 

3x200 

2 

3 

2 

7671 

7651 

f754/50 

3x200 

2 

3 

2 

7678 

7658 

f754/50 

3x200 

3 

4 

2 

7685 

7665 

f770/50 

3x200 

3 

4 

2 

7692 

7672 

f770/50 

3x200 

3 

4 

3 

7699 

7679 

f770/50 

3x200 

3 

4 

3 

7706 

7686 

f770/50 

3x200 

4 


3 

7713 

7693 

f770/50 

3x200 

4 


3 

7720 

7700 

f770/50 

3x200 

4 



7727 

7707 

f770/50 

3x200 

4 



7734 

7714 

f770/50 

3x200 


Main Survey 


A901a/F21 


A902/F22 



Figure 3. A cartoon illustrating the design of the observations for the main survey and the two cluster core fields (F21 and F22). 
Note that this illustration is valid in the case of a single galaxy, for which the wavelength is almost constant (for a given image there is 
wavelength variation across the field-of-view, see Eq. ^and Fig.^. A given galaxy is observed using a series of wavelength setups Xi and 
a 3-position dithered pattern. In each dithered position the wavelength changes slightly because the distance to the centre of the field 
changes. The order-sorting filter (OSF) separates the different tuned wavelengths. The observations are executed in ‘Observing Blocks’ 
(OBs). See http://www.nottingham.ac.uk/astronomy/omega/videos/ for videos with successive frames of a given galaxy. 


flat. The vignetted regions of the images are also trimmed 
at this stage. 

The final step carried out within OOPs is the sky sub¬ 
traction. Due to the wavelength variation across the FOV, 
sky OH emission lines appear as rings (see top panel of 
Fig.§. To correct for these features, each science exposure 
is artificially dithered on a 3 by 3 grid. The dither step size 
(a few arcseconds) is larger than the vast majority of the 
sources, but small compared to the ring features. Therefore, 


median-combining the dithered images leaves only the sky 
rings. These can then be subtracted from the original image. 
Fig. [^presents the first image of Field 21, with Ac = 7615 A. 
The upper panel shows the raw image while the bottom 
panel displays the final output from OOPs. 

As mentioned in Section l^we obtain three dithered ex¬ 
posures for each Ac (see FigTl^. This dithering introduces 
small shifts in the wavelength corresponding to each pixel. 
As we wish to retain this information, to improve our spec- 
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Z 


Figure 2. Spectroscopically determined redshift distribution of 
galaxies within 1.2 Mpc diameter apertures around the three 
main cluster cores, from a limited spectroscopic survey of the 
300 brightest galaxies in the region, taken with the AAT 2dF 
instrument. Arrows mark the positions of the brightest clus¬ 
ter galaxies. Dashed lines indicate the velocity range probed by 
OMEGA: the top histogram indicates the limits of the survey 
proper (fields 2, 3 and 5-20), while the ranges for the two pointed 
observations of the A901a and A902 cluster cores presented in 
this paper (fields 21 and 22 respectively) are specially tuned (cf. 
Table [l|. 



7600 

7575 

7550 

7525 

7500 

7475 

7450 

7425 

7400 


Figure 4. A mock reduced image with Ac = 7615 A, illustrating 
the radial wavelength dependence according to the colour bar 
(labelled by wavelength in A). The vertical line in the middle 
represents the 9.4 arcsec gap between CCDl and CCD2 and the 
axis are in pixels. 



Figure 5. Example of the first image of F21 with Ac = 7615 A. 
Top panel: raw image. Bottom panel: after OOPs processing. The 
red circle has a radius of 4 arcmin, within which the observations 
are free from contamination by other interference orders. 


trometry solutions are found using CCMAP and CCSETWCS 
with the COMBO-17 i7-band image as reference. In order 
to avoid unnecessary distortions, we work with the two sep¬ 
arate CCDs of OSIRIS throughout the analysis, obtaining 
astrometric solutions with rms accuracies ^0.1 arcsec. 


tral sampling, we do not coadd the dithered images. To cor¬ 
rect for cosmic rays, we therefore employ a similar routine 
to QZAP in IRAF applying a 3cr clipping with a fluxratio of 
0.1 (comparison value between the flux in the surrounding 
pixels and the flux in the candidate one to check whether 
the pixel has been hit by a cosmic ray). 

Bad pixels are corrected using the IRAF task FIXPIX. As- 


3.2 Wavelength Calibration 

The nature of the observations mean that the wavelength of 
a given object will depend on the central wavelength of each 
image as well as its distance to the optical centre. Both pa¬ 
rameters change from exposure to exposure. The variation of 
wavelength across the field-of-view for a given central wave¬ 
length is shown in Fig.|^ 

The wavelength calibration of the OSIRIS tuneable fil- 
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ter has been a question of debate (Mendez-Abreu et al. 


2011) because its accuracy is crucial for studies of extended 


sources, as well as for estimating redshifts of both point-like 
and extended sources. Fortunately, the sizes of our objects 
are small (the largest objects in F21 and F22 have a radius 
of ~ 7 arcsec), and hence the wavelength changes are neg¬ 
ligible across individual objects in the OMEGA sample (at 
most ~ 2A, i.e., much smaller than the spectral resolution, 
14AFWHM). The general form of the A calibration follows 
from Eq.[^ The optical centre is located on pixel Xq = 1051 
and To = 976 (CCDl) and Aq = -10 and To = 976 (CCD2) 
of the raw images, according to the OSIRIS website Q This 
corresponds to Xq = 772, Yb = 976 (CCDl) and Xq = —35 
and Yq = 976 (CCD2) on the images output from OOPs. 
Although there are other optical centres listed in the liter¬ 
ature ( Mendez-Abreu et al.[|2011 and the OSIRIS manual), 
the change they introduce is smaller than 0.2 per cent in 
wavelength. Using either of these values makes no difference 
for the purposes of the current paper. Here we are inter¬ 
ested primarily in the measurement of the Ka and [Nil] line 
fluxes, and not in the determination of accurate redshifts. 
However, in future papers, the combination of the images 
for all 20 OMEGA fields will allow us to evaluate the accu¬ 
racy of the wavelength calibration using a very large number 
of galaxies in the overlap regions (see Eig. 0 that are ob¬ 
served two or even three times. Eurthermore, the available 
AAT 2dE spectroscopic redshifts will provide an external 
cross-check. 


3.3 Flux calibration and zero-point estimation 


Different image sets (OBs) for a single field are not necessar¬ 
ily observed on the same night or with the same observing 
conditions. In order to perform consistent photometry, we 
therefore need to correct for differences in the seeing condi¬ 
tions, transparency and airmass. 

To match the seeing conditions our procedure is as fol¬ 
lows: we measure the PSE in all the images using unsat¬ 
urated stars with mi?,coMBO-i 7 < 19 (from the STAGES 
catalogue of Gray et al.|2009 ). We then convolve each image 
with a Gaussian function of the required width to produce 
images with the same PSE. Erom now on we use images 
matched to the worst seeing, which is 1.2 arcsec. These op¬ 
erations were performed using the IRAF tasks PSFMEASURE 
and GAUSSIAN. 

We perform photometric measurements using PHOT 
within IRAF on the stars in each individual OSIRIS image 
(typically, 15-25 stars per field) using an aperture of 2 arcsec 
(8 pixels) radius, sufficient to measure the total flux for these 
point sources. Because the wavelength at which each star is 
observed varies for different images, we obtain a wavelength- 
dependent zero-point zp{X) function. 


Zp(X) = m*W753f — OSIRIS , 


( 3 ) 


where m*w753f is the standard magnitude of that star, de¬ 
fined below, and m* osmis is the instrumental magnitude of 
the star measured in PHOT. 

To calibrate our fluxes we use the Vega magnitudes of 
the same stars in the closest GOMBO-17 filter, W753f, as 


^ http://www.gtc.iac.es/instruments/osiris/ 



A (A) 


Figure 6. Top panel: The wavelength-dependent zero-point, 
zp{X) = OSIRIS ~ 7n*w753f5 for 4 different dithered images 
in Field 21. Each set of ^^25 stars is shown with a different colour 
and symbol. We only show 4 of the total 42 images for clarity. In 
order to model the effect of the telluric line, we fit a polynomial. 
The best solution for each image set is shown by the lines. Bottom 
panel: The residuals with respect to the polynomial fits. 


standard magnitudes. This filter samples the range 7440- 
7620 A. Vega magnitudes were calculated from the photon 
fluxes tabulated in the GOMBO-17 catalogue. Our zero- 
point functions are shown in Eig.[^ Because there is a strong 
absorption feature caused by a double telluric sky line at 
A = 7633 A, the zero-point function has a strong wavelength 
dependency. 

In order to correct for the effect of the atmospheric tel¬ 
luric absorption, we fit this feature with a series of Gheby- 
shev polynomials. The shape of the absorption feature is 
constant, except for a scaling factor and offset that vary 
with the observing conditions. We therefore fit all the zero- 
point functions simultaneously, to obtain a general set of 
polynomial coefficients, together with an offset and scaling 
factor for each set of dithered images (14 in the case of Eield 
21, and 12 for Eield 22). We use a similar fitting technique to 
that described later, in Sec. |4.1| In Eig.j^we show the zero- 
points, as well as our fits and their residuals, as a function 
of wavelength. Eor clarity, we only show these for 4 images 
from the set of 14 dithered images. 

The scatter in the zero-point introduces additional 
wavelength-dependent errors that must be incorporated in 
our photometric uncertainties. We estimate these by tak¬ 
ing a 3cr-clipped standard deviation of the residuals in bins 
of 10 A. These zero-point errors are found to be always 
<0.1 mag. 


3.4 Galaxy photometry 

We perform aperture photometry on the galaxies with 
PHOT. We measure each galaxy using two apertures: a PSE- 
matched aperture of radius 1.2 arcsec (5 pixels, from now 
on referred to as Rrsf) and a total aperture (essentially 
the Kron radius; Bertin X Arnouts] |1996| ), with Rtot = 

2.5 X SMA. Here SMA is the second central moment of the 
light distribution, derived from AGS images with SExtrac- 
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TOR and contained within the STAGES catalogue (Gray et 


ah 2009). We use these two different aperture radii as ap¬ 


propriate throughout the paper, i^psp provides information 
about the nucleus of a given galaxy, while Rtot averages over 
its full extent. Fig. shows the spectra constructed for one 
of the galaxies in F21 using both i^psp or Rtot apertures. 

Since at different wavelengths the galaxies’ light distri¬ 
butions may change, we do not re-centre the aperture for 
each measurement, but instead fix the aperture position to 
the RA and Dec from the STAGES catalogue of [Gray et al.| 
(2009). The flux calibration is then applied using the zero- 
point zp{X) previously obtained via Eq. 


R^caiib(A) = mosiRis(A) — zp{X) . 


(4) 


When working in fluxes this becomes 

Ecaiib(A) = EbsiRis(A) X , (5) 

where iTaiib(A) is our calibrated flux in ergscm“^s“^A ^ 
and EbsiRis(A) is the flux measured with PHOT. 


3.5 Building the spectra 

After the flux and wavelength calibrations are performed, 
we are able to construct a spectrum for each of the galaxies 
in the observed fields. An example of such a spectrum can 
be found in the left panel of Fig. where we indicate the 
location of Ha and the two [Nil] lines with their respective 
rest-frame wavelengths. By design, the wavelength range of 
the spectra for cluster members should cover both Ha and 
[Nil] A6583, which will appear as emission lines in galaxies 
hosting star formation and/or AGN activity. In addition to 
finding spectra of galaxies with clear emission lines, we also 
expect to obtain spectra of galaxies with weak or absent 
emission (or absorption in Ha) lines. 


4 ANALYSIS 
4.1 Spectral fitting 


Our goal is to measure fluxes for Ha as well as the brightest 
[Nil] line at 6583 A. We also wish to obtain reliable con¬ 
straints for the continuum, in order to measure the equiv¬ 
alent widths {Wkcx and IT[nii]) of these lines. Due to the 
limited wavelength coverage and the moderate spectral res¬ 
olution and sampling of our spectra, this exercise becomes 
very challenging for weak emission-line galaxies (ELGs). 

To account for all the features present in our spectra, we 
fit them with a composite of three Gaussian functions and a 
linear continuum. Since the ratio of the fluxes of the two 
nitrogen lines ([Nil] A6583/[Nll] A6548) is fixed by atomic 
physics to a value of 3.06 (determined by the ratio between 
the transition probabilities of each line, [Osterbrock Fer^ 


land 1989), we only need one parameter to fit both their 


fluxes. We therefore fit, at most, six parameters: continuum 
intercept, continuum slope, redshift, flux of Ha, width of 
Ha, and flux of [Nil] A6583. 

In order to exploit all the information contained in 
our data, the fitting procedure must evaluate all possible 
combinations of the parameters in a thorough yet efficient 


way. Gonventional fitting techniques, such as Levenberg- 
Marquardt optimisation, can be sensitive to the initial val¬ 
ues, and can become trapped in local maxima in the param¬ 
eter likelihood space. These problems can be ameliorated 
by using Markov Ghain Monte Garlo (MGMG) techniques. 


In particular, we use the software emcee (Foreman-Mackey 
et al.||2013 ), which implements the Affine Invariant MGMG 


Ensemble sampler of Goodman & Weare (2010). 

MGMG methods aim to draw samples from the pos¬ 
terior probability density distribution, which describes how 
likely each possible set of model parameters is given the 
data. From Bayes Theorem, the posterior is given by the 
product of the prior distribution (which can be thought of 
as pre-existing constraints on the parameter values) and 
the likelihood function (which describes how closely the 
data correspond to the model, for a given set of parame¬ 
ters). MGMG algorithms explore the parameter space along 
a chain. At each step, new values for the parameters are 
proposed. These are accepted or rejected, depending on the 
probability of the proposed values relative to that of the 
current parameter set. If the probability is higher, then the 
proposed parameters will be adopted. However, if the proba¬ 
bility is lower, there remains some chance of the new values 
being accepted. In this way, the parameter space is thor¬ 
oughly explored, and ultimately sampled in proportion to 
the posterior probability density. 

The Affine Invariant MGMG Ensemble sampler uses 
multiple ‘walkers’ to produce a set of chains. These are 
started at different points of the parameter space. Each 
walker explores following a proposal distribution that de¬ 
pends on the current position of the rest of the walkers. In 
this way the parameter space can be mapped more efficiently 
than with the original Metropolis-Hastings algorithm. We 
also make use of parallel tempering, in which additional sets 
of walkers take significantly larger steps than the standard 
set. This eases communication between different maxima in 
the case of multimodal posterior distributions and generally 
speeds up convergence. 

The performance of the MGMG depends on two issues: 
the ‘initialization bias’ and the ‘autocorrelation in equilib¬ 
rium’ ( Sokal [1996 ). The first one refers to the need of reach¬ 
ing the equilibrium distribution and therefore ‘erasing’ any 
dependence on the initial conditions. This can be achieved 
by estimating the autocorrelation time in the burn-in phase, 
and producing chains long enough for them to move away 
from the initial conditions. The burn-in phase will then be 
discarded and not included in the posterior distributions. 
The second one implies that once the chains have reached 
equilibrium, we need to take enough samples to obtain a 
statistically-significant estimate of the errors. Therefore the 
performance of MGMG will strongly depend on the length 
of the chains and the number of walkers, but also on the 
election of the priors. To ensure a good performance for our 
analysis we use 300 walkers, a burn-in phase of 500 iterations 
(~10 times the autocorrelation time) which are discarded, 
and 500 iterations after the burn-in phase to sample the 
posterior distributions. 

An example of the performance of this algorithm using 
a four-parameter model (described below) for the spectrum 
of one of our galaxies is shown in Fig. Unless otherwise 
specified we summarise our fit results using the median value 
of the (marginalised) posterior distribution of each param- 
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Figure 7. Left panel: The flux-calibrated spectra for the central (fllled red circles) and total (open green circles) apertures of the bright 
(mji ~ 18.7) galaxy with STAGES ID 42713. The locations of Ha (Aq = 6562.8 A) and the two [Nil] lines (Aq = 6548.1 A and 6583.5 A) 
are indicated. Right panel: Image of galaxy 42713 at A = 7646 A, with the total and PSF apertures overlaid as solid green and dashed red 
lines, respectively. This image has not been convolved to the standard PSF although the fluxes are actually measured on the convolved 
images. The image is 40 pixels on a side, which corresponds to 10 arcsec or 28 kpc at 2 ; = 0.167. The colour bar indicates the flux in 
counts. 


eter, and estimate uncertainties from the 16th- and 84th- 
percentiles. 

A crucial element of this process is to set the priors 
that appropriately constrain the regions of parameter space 
to be explored. From what we already know of our model 
and data, we set the following priors: 


(i) ^[Nll] A6583 > 0; 

(ii) redshift consistent with the Ha line being within the 
wavelength range probed by our spectrum, with an addi¬ 
tional 20 A window on each side; and 

(hi) F[nii] A6583 < 3 (Tka + coutinuum). 


The first of these forces the flux of [Nil] A6583 (which, phys¬ 
ically, can never be negative) to be positive; the second as¬ 
sumes that the Ha line lies within the spectral range and 
the third sets an upper limit for the [Nil] A6583 flux such 
that it cannot be larger than 3 times the flux of Ha (which 
is the maximum value of [Nil] / Ha observed in the WHAN 
diagram Cid Fernandes et al.|2010| , including the possibility 
of up to 3 A equivalent-width of absorption in Ha. 

As mentioned above, we define our model using six dif¬ 
ferent parameters. This number can be reduced if we assume 
some properties of our spectra. First, one can assume that 
the slope of the continuum is roughly the same for all the 
galaxies, which is reasonable considering the short wave¬ 
length range covered by the OMEGA spectra (^ 200 A). 
Secondly, the width of the emission lines will be limited by 
the instrumental resolution and only substantially broad¬ 
ened in cases of strong AGN (1000 - 25000 kms“^). Due to 
the low spectral resolution and the faintness of some of the 
objects, adding these constraints will improve the reliability 
of our estimates. 

In order to find a fixed value for the slope of the contin¬ 


uum and for the instrumental width of the lines, we first per¬ 
formed six-parameter MGMG runs considering only those 
galaxies that show clear emission lines. We then select galax¬ 
ies with less than 15 per cent uncertainty in the estimation 
of the flux and width of the Ha line, and the median value of 
the redshift consistent with Ha and [Nil] A6583 being within 
the wavelength range. In total we find 28 galaxies that ful¬ 
fill these requirements. For these galaxies. Fig. shows the 
distribution of the line-width and the slope of the contin¬ 
uum. The continuum slope is consistent with zero, while the 
distribution of the width of the lines shows a clear peak at 
7.5 A. These findings allow us to confidently fix these two 
parameters and use only four-parameter models to fit the 
spectra. 


4.2 Defining the OMEGA sample of ELGs 


After fitting the Ha and [Nil] lines as described above, in this 
section we describe the different criteria we use to define the 
sample of emission-line galaxies we will use in the remain¬ 
der of the paper. We start with all objects in the STAGES 
catalogue with 16 ^ tur ^23.0 that are detected in both 
GOMBO-17 (combo_flag > 1) and in HST (stages_flag > 
1). At this stage we do not require them to also be flagged 
as clus ter members according to the definition of [Gray et al.| 
( |2009| lp] 


^ Membership of the A901/2 cluster structure is determined tak¬ 
ing into account the COMBO-17 photometric redshifts and their 
uncertainties. Given the high accuracy of these redshifts, the pro¬ 
cedure maximises completeness without introducing unacceptable 
contamination. Details ara given in Sect. 3.3 of |Gray et al.| ( [^09| >. 
See, in particular, their Fig. 13 and Equation 8. 
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Figure 8. Parameter-parameter likelihood distributions with one- (solid), two- (dashed) and three- (dash-dotted) sigma contours obtained 
for the four-parameter model of the Rtot aperture spectrum of galaxy 42713 (the same galaxy shown in Fig.[^. On the top right we 
show the observed spectrum (blue points) together with the model fits (red lines). We plot only 1/100 of all the individual MCMC fits 
for clarity, but they are all very similar. The median values of the parameters’ probability distribution are plotted as a white cross on 
the two-dimensional distributions. The parameters shown are the continuum flux density, the redshift, and the integrated fluxes of the 
Ho; and [Nil] A6583 emission lines. Fluxes are given in c.g.s units. 


We use measurements made within the total aperture, 
i7tot, as defined in the previous section. Our final sample 
of ELGs in F21 and F22 contains 124 objects fulfilling the 
following selection criteria: 

(i) relative error(WHa) < 0.15; 

(ii) probability of detecting the Ha line in emission, 
P{Fiia>0) > 99.7; and 

(iii) the Ha line falls within the wavelength range probed 
by our observations. 


Out of these 124 objects, 103 have been morphologically 
classified by the STAGES collaboration (see, e.g., Wolf et 
al.|2009 Malt by et al.|2010 ). The galaxies without morpho¬ 


logical classification have Zphot >0.4, the STAGES cut-off 


for visual classification. Of the 124 ELGs, 72 are classified as 
clusters members in the STAGES catalogue based on their 
GOMBO-17 photometry (combo_flag > 3). 

Eor the parts of our analysis that make use of both Ha 
and [Nil] measurements (e.g.. Sect. |5.1| , we define a sub¬ 
sample requiring both Ha and [Nil] A6583 to be within our 
wavelength range. This subsample contains 82 ELGs with 
both i^tot and i^psF aperture measurements. 


4.3 Detection efficiency and limits 

Our ability to detect Ha emission depends on a complex 
interplay between the brightness of a galaxy, its size, and 
its Ha flux and equivalent width. Eor instance, for a given 
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Figure 9. Distribution of the width (a) of the spectral lines and 
the continuum slope for a sample of 28 galaxies with very clear 
emission lines. As a reference, the continuum level of a mji = 20 
galaxy is ~ 5 x 10“^® ergs cm“^ 


Figure 10. The distribution of the Ha flux for the 124 ELGs in 
our sample, where the error bars are 1/y/N. 


Ha flux, galaxies with stronger continuum will be harder 
to detect. Figs, and m show the distribution of the Ha 
fluxes and equivalent widths (Wna) for the 124 ELGs in 
our sample. We are able to detect galaxies with Ha fluxes 
as low as ~ 1.3 x 10“^^ ergcm“^s“^. Given the relatively 
small sample presented in here (only galaxies in F21 and 
F22), we leave for a later paper a detailed discussion of the 
Ha flux and equivalent width distributions and their de¬ 
pendence on the intrinsic properties of the galaxies and the 
environment. Reassuringly, the Hfoa distribution we obtain 
is broa dly compatible with the one found by |Bafogh et al.| 
(2004). For equivalent widths larger than ~ 20-30A the dis- 
tribution declines steeply, following roughly a power law. At 
lower equivalent widths the distribution flattens out until 
incompleteness probably kicks in below ~ 5A. When inter¬ 
preting our data it is important to bear in mind that the 
sensitivity of our survey declines for lower equivalent widths. 

Due to the complex nature of the Ha detection process, 
it would be very difficult to model accurately the detec¬ 
tion limits and the completeness of our sample from first 
principles, instead, we follow an empirical approach that 
takes advantage of the wealth of data available from the 
STAGES database. Ha emission is expected to be present 
in blue, star-forming galaxies. Moreover, many red spirals 
have been shown to have a certain degree of star forma¬ 
tion ( Wolf et al.|2005| ). Here we compare our Ha detections 
with the numbers of blue-cloud and red spiral galaxies (star 
forming candidates) in the STAGES catalogue in order to 
evaluate how efhcient the OMEGA survey is at detecting 
star-forming galaxies and to assess the sample completeness. 

fn the top panel of Eig.^jwe show the rriR distribution 
for the galaxies in our sample where Ha has been detected 
(green). Eor comparison, we overplot in yellow the distribu¬ 
tions for all GOMBO-17 cluster galaxies in the same field- 
of-view and in blue those identified as star forming (SE). SE 
galaxies are either blue-cloud or red-spiral galaxies ( [Gray et 


^|2009 ). The dashed histograms show the same GOMBO- 


17 samples corrected for completeness and contamination 



logCWnJ 

Figure 11. The distribution of Wna for fhe 124 ELGs in our 
sample, where the error bars are I/VTV. For equivalent widths 
larger than ~ 20-30A the distribution declines steeply with in¬ 
completeness kicking in below ~ SA. We have fitted a line with 
slope -1.5 and an intercept of 3.8. 


(see below). One caveat here is that we do not separate the 
so-called “retired” galaxies ( Gid Fernandes et al.|20lT Yan 


Blanton] 2012), which also have Ha emission, from the 
star-forming ones. More details are given in Sect. |5.l1 when 
we discuss AGN activity. 

The parent sample of GOMBO-17 galaxies in F21 and 
F22 contains ~ 1400 objects, with 326 galaxies expected 
to be cluster members (listed in the STAGES catalogue as 
having 16.0 ^ ttir ^ 23 and combo_flag > 3) with stellar 
masses 9.0 ^ log (M*/Mq) ^ 11.5. Of these, 203 are either 
red spirals (sed_flag_cl = 2) or belong to the blue cloud 
(sed_flag_cl = 3), both thus considered to be star-forming 
galaxies, based on their GOMBO-17 SEDs. We securely de¬ 
tect Ha in 124 objects. Erom this we estimate that we are 
able to detect Ha in ~ 40% of the cluster members and 
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in ~ 60% of the cluster star-forming galaxies more massive 
than 


As described in Gray et ah (2009), the cluster mem¬ 
bership was determined using COMBO-17 photometric red- 
shifts. Despite the relatively high accuracy of these redshifts, 
the cluster sample suffers from a certain degree of incom¬ 
pleteness and contamination, particularly at faint magni¬ 
tudes. Therefore, the “true” (corrected for incompleteness 
and contamination) number of cluster galaxies A^cmcir can 
be estimated as 




A7;mcl7(l COntfrac) 


(6) 


cornpfrac 

where iVcmcir is the observed number of cluster mem¬ 
bers based on their photometric redshifts. In this equa¬ 
tion, and cont^rac are the magnitude-dependent 

completeness and contamination fractions respectively (see 
Fig. 14 of Gray et al.|20(39 ). Similarly, the “true” number of 
SF cluster galaxies iVsFcmcir is 

■^jc _ A^SFcmcl7(l COTT-tfrac) 

^''SFcmcl7 — -5 VG 

compfrac 

where iVsFcmci 7 is the observed number of SF cluster galax¬ 
ies. 

The bottom panel of Fig. shows as green stars the 
fraction of cluster galaxies detected in Ha, defined as 

det 


/na det — 




(8) 


and as blue circles the fraction of SF cluster galaxies de¬ 
tected in Ha 


/sFHadet = 


A^Hac 


(9) 


We effectively detect Ha for all SF cluster members 
brighter than ttir ~ 20. Fainter than this, the complete¬ 
ness drops off, as expected - only galaxies that are bright in 
Ha are detected. Note that the apparent drop for galaxies 
brighter than rrir ~ 19 is due to low-number statistics in 
the brightest bins. As a fraction of all cluster-members, the 
completeness does not depend strongly on ttir since at lower 
masses a greater fraction of cluster galaxies are star forming. 
Note that /det > 1 means that we are detecting Ha in more 
galaxies than are expected to be star forming (based on the 
sed_flag_cl flag). This reflects both the uncertainty in the 
GOMBO-17 SED classification and the high sensitivity of 
our survey. 


5 RESULTS 


As a preview of the kind of scientific issues that we will be 
able to address using OMEGA data, in this section we iden¬ 
tify star-forming galaxies and AGN from our ELG sample, 
and investigate their properties, including mass and mor¬ 
phology. This is done by cross-matching our OMEGA obser¬ 
vations with the STAGES catalogue ( Gray et al.||2QQ9| . In 
addition, we present the Ha luminosity function (LE) in the 
densest regions in order to understand how comparable our 
Ha SERs are to literature values of other sources. We finally 
compare our Ha-derived star-format ion rates (SERs) to UV 
and IR ones from Gray et al. (2009). In subsequent papers 


we will present a full analysis for the complete dataset. 



Figure 12. Top panel: the uir distribution of galaxies deter¬ 
mined to be cluster members using the COMBO-17 photometric 
redshifts (sample denoted cmcl7, in yellow), and the subset of 
these that are expected to be SF (galaxies classified as blue-cloud 
or red spiral, in blue). The dashed histograms indicate the same 
cmcl7 and SF cmcl7 samples, but corrected for completeness and 
contamination as discussed in the text. The green square-filled 
histogram is the distribution of galaxies for which Ha has been 
detected in our OSIRIS spectra. Bottom panel: The Ha detection 
rate, expressed as a fraction of the full sample of cluster members 
(/Ha det; green stars) and as the fraction of SF cluster members 
(/sFHadet; blue circlcs), plotted as a function of ruR. See text 
and equations and for details. 


5.1 AGN and star-forming galaxies 


In Fig. 


13 


_ we show the [Nil] / Ha vs. Wna (WH AN; pd| 

Fernandes et al.|[2010 |Gid Fernandes et al.||20ll] ) diagram 
for our ELGs. Here we use the Rpsf aperture, as this pro¬ 
vides the cleanest measurement of the nuclear emission of 
the galaxy. The WHAN diagram is an alternative economical 


form of the widely-used BPT diagnostic diagram (Baldwin 
et al.||1981 ). Gid Fernandes et al. (2010) and Gid Fernan¬ 
des et al. (2011) have convincingly shown that the ratio of 


[Nil] to Ha is truly powerful in separating AGN from SF 
galaxies, especially when the equivalent width of Ha is also 
considered. 

Using the probability density distributions obtained for 
the continuum level and fluxes in Ha and [Nil], we can ob¬ 
tain the probability density distributions of the Wna and the 
[Nil] / Ha flux-ratio. Assuming standard values to differenti¬ 
ate sources of ionization in the WHAN diagram (cf. Fig. 13), 
we can therefore calculate the probability for each galaxy 
in our sample being dominated by star-format ion or AGN 
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activity. Following this, 10 objects have a greater than 99.7 
per cent probability 3a confidence) of hosting a dominant 
AGN (Seyfert; red circles), while 16 have a similar proba¬ 
bility of being SF-dominated (blue stars). The remaining 56 
galaxies (green points) have less certain classifications due 
to the observational uncertainties on [Nil] / Ha. Obviously, 
it is possible to classify galaxies into SF and AGN with less 
certainty by lowering the probability threshold. This would 
increase the size of the SF and AGN samples, but the con¬ 
fidence for each individual object would be reduced. Nev¬ 
ertheless, since we have robust probability distributions for 
the line ratios, the statistical properties of such samples can 
still be analysed rigorously. It is important to bear in mind, 
however, that in reality there is no hard division between 
SF- and AGN-dominated systems. Better observational er¬ 
rors would place the galaxies more accurately on the WHAN 
diagram, but since the boundaries are somewhat fuzzy (and 
depend on other physical parameters such as metallicity and 
inter-stellar medium density), the classifications would, in 
many cases, remain uncertain: a significant fraction of the 
green points may well contain mixed sources of ionisation. 


One interesting population revealed by the WHAN di¬ 


agram is the so-called “retired” galaxy population (Gid Fer- 
nandes et al.|[20TT] |Yan fc Blanton] 12012) [Stasihska et al. 


2015| . They are located in the bottom-right part of such di¬ 


agram, which is also occupied by LINER-type AGN. These 
retired galaxies have stopped forming stars, their gas is be¬ 
lieved to be ionized by hot low-mass evolved stars, and they 
do not have a detectable contribution from an AGN. Fur¬ 


thermore, their Ha emission is often extended (Belfiore et al. 


2014). This population will be addressed in detail in a sub¬ 


sequent manuscript, when we analyse the resolved emission¬ 
line images. 


In the following we explore how the properties of the 
SF galaxies and AGN hosts that we have identified com¬ 
pare to one another and to the overall population of ELGs. 
We consider here their total ttir, morphologies and stellar 
masses. 


In Eigs.p^andp^we present histograms of ttir and stel¬ 
lar masses for the ELGs with both Ha and [Nil] A6583 secure 
detections. The distributions of the secure AGN (Seyferts) 
and SE objects are also shown. Galaxies hosting AGN are 
brighter in the R-band and more massive than SE galaxies. 
Moreover, the overall behaviour of the Ha-detected galax¬ 
ies appears to mimic the mass and ttir distribution of SE 
galaxies, suggesting that the majority of the the galaxies 
with uncertain classifications are SE-dominated. 

In Eig. we show the distribution of morphologies. 
Galaxies determined to be hosting AGN with ^ 3a confi¬ 
dence are typically of earlier-type than those for which the 
dominant source of ionisation is SE or uncertain. 


In Eig. we show the distribution on the sky of the 
galaxies in the E21 and E22 regions down to ttir < 23 that 
are classified as cluster members in the STAGES catalogue. 
Galaxies with more than 3a conhdence of being AGN- or SE- 
dominated (see Eig. 13) are denoted by red and blue filled 
circles respectively, whereas other Ha secure detections are 
indicated by filled green circles. The centres of A901a and 
A902 are marked with black crosses and are defined as the 
location of the brightest cluster galaxies. We note that secure 
SE and AGN galaxies are located towards the outer parts of 



Figure 13. A diagnostic plot of [Nil] / Ha vs. Wna-, the so-called 
WHAN diagram ( |Cid Fernandes et al.|[?Q10| ). The vertical and 
horizontal lines separate regions with different sources of ioni¬ 
sation ( |Cid Fernandes et al.||2Q10| |2011| ). We only show detec¬ 
tions using RpsF aperture measurements where both Ha and 
[Nil] fall in the OMEGA probed wavelength range: a total of 
82 ELGs. From these we find 10 secure (^ 3f7 confidence) AGN 
hosts (Seyferts; red circles) and 16 star-forming galaxies (SF; blue 
stars). The green points represent galaxies for which the classifi¬ 
cation is uncertain and/or contain mixed sources of ionisation. 



Figure 14. The distribution of the R-band magnitudes from 
COMBO-17 for ELGs with detected [Nil] A6583 and Ha (green 
line), as well as for the subsamples of high-confidence AGN 
(Seyferts; red histogram) and SF objects (in blue). See text for 
details. 


the clusters and are completely absent in the central regions, 
especially when looking at F22. 


Gilmour et al. (2007) presented a catalogue of XMM- 


detected point sources in the A901/902 field (see their ta¬ 
ble 12). From the 12 X-ray sources that they identify as 
likely AGN within the cluster, 4 objects fall in the FOV 
of F21 and F22. Only one of them is part of our ELG 
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Figure 15. Stellar mass distribution for ELGs with detected 
[Nil] A6583 and Ha, together with the AGN and SF subsamples. 
Object selection and lines as in Fig. |14| 



Figure 16. The distribution of morphologies for FLGs with de¬ 
tected [Nil] A6583 and Ha as well as the subsamples of the high- 
confidence 3(7) AGN and SF objects. 


sample with securely-detected Ha (ID 11827). However, this 
galaxy is not classed as an AGN because the [Nil] A6583 
line used in the classification process falls outside the wave¬ 
length range probed by the F22 observations. Reassuringly, 
the [Nil] A6548 line is detected and, based on its intensity, 
this object has a 93.5% probability of being an AGN. More¬ 
over, preliminary analysis of an overlapping field (which 
samples a slightly different wavelength range) indicates that 
[Nil] A6583 is clearly detected for this object, indicating that 
it is probably an AGN. For a second X-ray AGN (ID41435) 
neither Ha nor [Nil] fall in the wavelength range probed 
by our observations. For the remaining two X-ray AGNs 
(ID 12953 and ID44351) our observations sample both Ha 
and [Nil] but their Ha lines are below our detection prob¬ 
ability threshold (P(FHa > 0) < 99.7; see section 4.2) and 
therefore do not make it into our ELG sample. In summary. 


one out of the four X-ray AGN that we could have detected 
is also highly likely to be an optical AGN. The others have 
either too faint emission lines or do not have Ha or [Nil] 
in our probed wavelength range. It is possible that these 
objects may be heavily obscured Gompton-thick AGN and 
therefore would not be expected to show optical signs of nu¬ 
clear activity. On the other hand, the fact that none of the 
10 highly likely optical AGN found in the OMEGA survey 
are associated with a Gilmour et al. (2007) X-ray source 
clearly shows that OMEGA is able to detect a large number 
of optical AGN with relatively weak or no X-ray emission. 


5.2 The Ha Luminosity Function 

The Ha luminosity function of star forming galaxies is an im¬ 
portant measurement for comparison among different stud¬ 
ies and environmental dependencies. The luminosity of the 
Ha nebular emission line is a SFR indicator with a direct 
physical connection to short-lived, massive stars. In order 
to calculate the luminosity function of Ha, line fluxes are 
converted to luminosities via the following: 

LyLcx — 477 Dl ihia , (10) 


where Dl is the luminosity distance (795 Mpc at z = 0.167). 

In Fig. we show the Ha cumulative luminosity func¬ 
tion (LF) for the highest density regions of A901/2, using 
our Ha detections from F21 and F22 as described in the 
previous sections. The top axis shows the SFR, assuming 
1 mag average extinction, according to: 


SFRh« [Mgyr-^] = 2 x [ergs”^] 


( 11 ) 


following |Kennicutt| (19981. 

We use the results of Sect. 14.31 in order to correct the 
OMEGA Ha LF for incompleteness. We parametrise the 
dependence of /sFHct det on the R-band magnitude ruR as 


/sFHa det 


—0.24mi^ + 5.7 if ttir ^ 20 

1.0 if tur < 20 


( 12 ) 


We consider our sample to be essentially complete for ttir < 
20 (see bottom panel of Figp^. To correct for incomplete¬ 
ness, we weight the OMEGA Ha LF by l//sFHadet- 

For comparison we show the A1689 luminosity func¬ 
tion (at z = 0.18) from Balogh et al. (|2002| as w ell as the 
field Ha LF from GAMA ( Gunawardhana et al.|[20T3 ). For 
the field Ha LF we use the best-fitting parameters of the 


Saunders function (Saunders et al. 1990) that Gunaward- 
|hana et al. (2013) obtain for the redshift bin z = 0.1-0.2 
(logL* = 34.55, logC = -2.67, a = -1.35, a = 0.47) to 
reconstruct the corresponding normalised cumulative frac¬ 
tion. For A1689 we use the Schechter functio n parameters: 
a = —0.1 and L* = 10^° ergs“^ given by Balogh et al. 


(2002), to reconstruct the cumulative Ha LF. Gompared to 
the field sample from [Gunawardhana et aT (2013), the dens¬ 
est parts of A901/2 contain galaxies with lower, on average, 
SFRs. Gonversely, A901/2 hosts brighter Ha emitters than 
A1689. This is not a surprise, as A1689 is a more evolved and 
denser cluster than the densest regions of A901/2 ( jLemze ^ 


aL]|2009 Alamo-Martmez et al.||2013 ). We conclude that 


the densest regions of the A901/2 system star formation has 
been significantly suppressed compared to the field, but not 
as much as in denser, more evolved, cluster environments. 
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Figure 17. The distribution on the sky of the cluster galaxies (as determined by COMBO-17; denoted cmcl7 in the legend) in the 
F21 and F22 regions down to mR < 23. ELGs (where Ho; has been securely detected) are represented by green filled circles. Galaxies 
classified as AGN or star-forming with ^ 3a confidence are shown as solid red circles and blue stars respectively. The centres of A901a 
and A902 (determined from the brightest cluster galaxies) are marked with black crosses. The large dashed circles show the r2oo radii of 
A901 and A902 derived from the weak leasing analysis of |Heymans et aT] ( |2008| >. 


5.3 Star Formation Rates 

The SFR-stellar mass relation is a useful diagnostic for the 
star-format ion properties and history of a galaxy population. 
In this parameter space there is a reasonably well-defined 
main sequence whose properties have been characterised at 


different redshifts (see e.g. Whitaker et al. 2012). Several 


tracers across the electromagnetic spectrum may be used 
to estimate the star formation rate of a given galaxy. The 
STAGES catalogue of [Gray et al. (2009) provides UV and IR 
SFRs for the A901/2 region. Rest-frame UV light originates 
mainly from massive stars and thus directly traces young 
stellar populations (with main sequence lifetimes < 10® yrs); 
however, it is strongly attenuated by dust. The SFR from 
the UV continuum is often dubbed ‘current SFR’. The in¬ 
terstellar medium associated with a star-forming region can 
be quite dusty, so a significant fraction of the UV lumi¬ 
nosity produced by massive young stars can be absorbed. 
This radiation heats the dust and is re-emitted in the far 
infrared. Moreover, the interstellar medium around young 
massive stars is ionised by the Lyman-continuum photons 
they emit. The recombination of this ionised gas produces 
Ha emission-lines. The SFR from the Ha nebular emission is 
often dubbed ‘instantaneous SFR’, tracing stars with main 
sequence lifetimes of < 2 x lO^yrs (Mo, van den Bosch and 


White|2(HT ). In the following figures, we adopt an Ha SFR 
with 1 mag of extinction (see Eq. 11). At this point, we do 
not attempt to subtract the AGN contribution to the SFR 
from the Ha fluxes. Resolved star formation properties will 
be studied in a subsequent manuscript. 

We compare our Ha SFR against those from other avail¬ 
able tracers in Fig. The UV SFRs are based on the UV 
luminosity estimated using the 2800A rest-frame luminosity 
from GOMBO-17, which has a complete coverage around the 
A901/2 multi-cluster system. The UV SFRs have been cor¬ 
rected for extinction as part of the SED-fitting procedure. 
Details can be found in Gallazzi et al. ( 20Q9| and [Wolf et| 
al. (2009). The IR SFRs are determined from Spitzer/MIPS 
24qm data. In the cases where there was no detection in 
the IR, upper limits to the SFRs were estimated ( [Wolf et ah] 
|2009| ). The green datapoints (denoted ‘UV+IR’ in the figure 
legend) represent the composite SFR measurements. Note 
that these points, which account for the obscured and unob¬ 
scured star formation, are, on average, ~ 0.3 dex above the 
one-to-one line. This difference may be due to the Ha SFR 
estimates missing heavily-obscured SF, but it could also be 
indicative of a declining SFR, given the different timescales 
of the indicators. In the case of the SFRs from the UV alone, 
the correlation with the Ha SFR is stronger, which is not 
surprising as the UV and Ha are sensitive to comparable 
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Figure 18. The cumulative Ho: luminosity function for the high¬ 
est density regions of A901/2 corrected (solid green line) and un¬ 
corrected (solid yellow line) for contamination and completeness. 
For comparison, we show the Ha LF from the cluster A1689 at 
2 ; = 0.18 (blue dash-dot line; |Balogh et al.|20Q2| ) and the field from 
the GAMA survey at 2 : = 0.1-0.2 (red dashed line; | Gunaward- 1 
|hana et al.|2013| ). The densest regions studied by OMEGA have 
SFRs that fall in between those of the cluster A1689 and the field. 



Figure 19. The logSFRHa determined from Rtot aperture mea¬ 
surements compared with other estimates of log SFR. See legend 
and text for details on the meaning of the different datapoints. 


SF timescales, with the UV tracing slightly longer ones than 
Ha. However, there is also a clear ~ 0.3 dex offset since the 
effect of extinction is significantly stronger in the UV than 
in the red. 

In Fig. we show the stellar mass vs. the Ha- 
determined SFR. Galaxies reliably classified as AGN and 
SF objects (cf. section 5.1 and Fig. 13) are indicated by 
red circles and blue stars respectively. At a given stellar 
mass, only the galaxies that are brightest in Ha can be reli¬ 


ably classified as SF or AGN, while the rest (green squares) 
have either mixed ionisation sources or uncertain classifica¬ 
tion, but they are expected to be dominated by star forma¬ 
tion (section 5.1). For comparison, the solid line shows the 
parametrized SFR-mass relation derived from equations 1, 2 
and 3 of |Whitaker et al.| ( |2012 ) for star-forming field galaxies 
at the redshift of A901/2, 


log SFR = 0.68 logM* - 6.56. 


(13) 


The stellar masses used by Whitaker et al. are de¬ 
rived from medium-band photometr y using|Bruzual Ghar-] 
lot (2003) models that assume a Ghabrier (2003) initial 


mass function (IMF), solar metallicity, and the Galzetti et 


al. (2000) extinction law. The stellar masses that we use 
were estimated by |Borch et al. (2006) using the GOMBO- 
17 photometry in conjunction with a template library de¬ 


rived using the PEGASE stellar population models (Fioc 


& Rocca-Volmerange 1997) with solar metallicity and a 


Kroupa (2001) IMF. The different IMFs used by Whitaker 


et al. and Borch et al. would only change the stellar masses 
by ~ 10%. Moreover, the stellar masses derived in these 
works are quantitatively consistent with those derived us¬ 


ing a simple colour-based stellar mass-to-light ratio (Bell et 
al.|2003 ). We are therefore confident that the stellar masses 


used here can be directly compared, and that any systematic 
differences will not affect our conclusions. The SFRs esti¬ 


mated by Whitaker et al. (2012) are based on UV+IR indi¬ 
cators, and are therefore expected to be on average ~ 0.3 dex 
higher, than the Ha-based ones. This has to be taken into 
account when comparing our data with that of Whitaker 
et al. Nevertheless, even if we consider such offset and any 
likely systematic uncertainties, it is clear that most of our 
Ha-detected galaxies would still fall clearly below the field 
SFR-mass relation, consistent with the differences found be¬ 
tween our Ha LF and the field one (Fig. p^. This confirms 
that many of the star-forming galaxies located in the centres 
of the A901/2 clusters appear to have their SFR significantly 
suppressed relative to the field. This suppression seems to 
occur at all stellar masses. 

At first sight, this seems to contradict previous results 
which suggest that the average SFR in SF galaxies remains 
roughly constant with environment, and it is only the frac¬ 
tion of SF galaxies that changes (Balogh et al. 2004| Ver- 


dugo et al.|2Q08 Poggianti et al.pOQS Bamford et al.pOQS ). 
However, this contradiction is only apparent. The relative 
fraction of SF and non-SF galaxies depends on the depth of 
the survey. Because we are able to detect very low Ha fluxes 
(and thus very low levels of star formation), our sample con¬ 
tains Ha-emitting galaxies that would have been missed by 
most other surveys. Therefore, strongly suppressed galaxies 
still appear as Ha emitters in our sample. Gluster galaxies 
with strong star formation (exemplified by the blue stars in 
Fig. 20) have levels of SF comparable with those found in the 


field. If only galaxies with strong star formation are detected 
in a survey, it would appear that their SF has not been sup¬ 
pressed. Since galaxies with very low levels of star formation 
would not have been detected in Ha in most other surveys, 
they would have been (erroneously) classified as non-SF and 
one would conclude that only the relative fraction of SF and 
non-SF galaxies has been affected by the environment. 


Wolf et al. (2009) reached the same conclusion based 


on restframe UV and 24/i-data of Abell 901/2. They found 
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Figure 20. Stellar mass vs. SFR(Hq;) determined from Rtot aper¬ 
ture measurements. The galaxies with a Scr confidence of being 
AGN- or SF-dominated, with selection based on the PSF aper¬ 
ture (see Fig |13| > are denoted by red circles and blue stars re¬ 
spectively. The dashed-line is a SFR-mass relation for field star¬ 
forming galaxies at z = 0.167 ( [Whitaker et al.|2Q12| >. Comparing 
our secure Ha detections to the dashed-line representing the field 
galaxies indicates a suppression of star formation in these high 
density regions. 


a large population of lower-SFR galaxies, mostly red spiral 
galaxies seemingly in transition from normally star-forming 
galaxies to quiescent ones. This population had higher FIR 
luminosities than expected from their broad-band optical 
colour alone, suggesting that the specific SFR of the aver¬ 
age red spiral at a given stellar mass was only three times 
lower than the SSFR of the average blue spiral at similar 
mass. An earlier study by Wolf et al. (2005) had already re¬ 
vealed this population from optical photometry alone, using 
mediumband imaging to discriminate between red colours 
due to high age in the absence of dust and those due to 
moderate age in the presence of moderate amounts of dust, 
thus differentiating between what they called ’old red’ galax¬ 
ies and ’dusty red’ galaxies. The average specific SFR esti¬ 
mated from [On] emission in their ’dusty red’ population 
was four times lower than that in blue galaxies. While this 
paper investigates only a small part of the same cluster, we 
find a similar transition population in SFRs estimated from 
Ha now. We are thus suggesting that this cluster contains 
a rich population of moderately star-forming galaxies fill¬ 
ing the previously considered gap between non-star-forming 
and star-forming galaxies. This view is now supported by 
evidence from a range of tracers, from UV and 24/i, from 
[Oil] and Ha-emission, and from a detailed analysis of the 
optical SED. A detailed study of the environmental depen¬ 
dence of the SF for the full galaxy sample will be presented 
in a subsequent paper. 


6 SUMMARY AND CONCLUSIONS 

In this paper we have introduced the OMEGA (Osiris Map¬ 
ping of Emission-line Galaxies in A901/2) survey. The ul¬ 
timate goal of this project is to study star formation and 


AGN activity across a broad range of environments at a sin¬ 
gle redshift. Using the tuneable-filter mode of the OSIRIS 
instrument on the 10.4m GTG we have targeted Ha and 
[Nil] emission lines over a ~ 0.5 x 0.5 deg^ region covering 
the z ~ 0.167 multi-cluster system containing Abell 901 and 
902. 

Based on observations of 2 out of the 20 fields covered 
by this survey, this paper outlines the data reduction proce¬ 
dures followed and describes the analysis we have developed 
in order to reliably measure the fluxes and equivalent widths 
of the Ha and [Nil] lines. We present the quantitative crite¬ 
ria used to define a robust sample of emission-line galaxies 
and analyse its properties in terms of detection efficiency, 
detection limits, and completeness. This paper also presents 
some examples of the kind of science that can be done with 
the survey. The scientific results we discuss are based on 
only 10% of the surveyed area, corresponding to the highest 
density regions of the A901/2 structure. An analysis of the 
full survey is left for future papers, but interesting results 
are already emerging. 


• In these two fields we detect 124 ELGs down to Ha 
fluxes ~ 1.3 X 10“^^ ergcm“^s“^, corresponding to Ha lu¬ 
minosities ~ 10^^ ergs“^ at the distance of the A901/2 sys¬ 
tem. This allows us to probe integrated SFRs as low as 
~ 0.02 Mq yr“^. We estimate that we are able to detect 
Ha emission in ~ 40% of the cluster members with stellar 
masses above ~ IO^Mq, and ~ 60% of the galaxies expected 
to be forming stars according to their SEDs in the same mass 
range. 

Using the WHAN diagram ( |Gid Fernandes et al.||2010 


2011), which plots Wna vs. the [Nil]/Ha ratio, we are able to 


identify 10 bona-fide AGN in the Abell 901 and 902 clusters. 
The AGN hosts tend to be brighter and more massive than 
the star-forming galaxies. Moreover, the AGN host galax¬ 
ies usually have earlier morphologies than the general star¬ 
forming galaxy population. 

• The emission-line galaxies we detect, both star-forming 
and AGN, tend to avoid the central highest-density regions 
of the clusters. 

• The Ha luminosity function of the galaxies in the dens¬ 
est regions of Abell 901/2 contains brighter Ha emitters 
than that of more massive and richer clusters at comparable 
redshifts (e.g., Abell 1689 at z = 0.18). However, our sample 
contains a lower proportion of very luminous Ha emitters 
than the field galaxy population at similar redshifts. This 
indicates a progressive suppression of the star-formation ac¬ 
tivity in galaxies when we move from the field to clusters of 
increasing richness/mass. 

• The SFR-stellar mass relation we observe in the dens¬ 
est regions of A901/2 falls well below the relation for field 
galaxies found by Whitaker et al. (2012) at similar redshifts. 
In other words, the average specific SFR (or SFR per unit 
stellar mass) is significantly lower in these clusters that in 
the field, confirming the star-formation suppressing effect of 
the clusters. 

• The SFRs determined from the Ha luminosities corre¬ 
late very well with the SFRs estimated from the rest-frame 
UV photometry, but are typically a factor of ~ 2 higher. 
The strong correlation indicates that the UV light and Ha 
trace roughly comparable star-format ion timescales but the 
effect of extinction is significantly stronger in the UV than 
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in the red. Ha-determined SFRs also correlate well with 
those determined combining the UV and IR emission, but 
the offset is now reversed: the Ha SFR estimates would miss 
heavily-obscured star formation and are thus a factor ~ 2 
smaller. Some of the difference could also be due to a declin¬ 
ing SFR, since Ha and the IR probe different star-formation 
timescales. 

We have found that the OSIRIS tuneable hlter is a 
highly effective instrument for mapping the emission lines 
of galaxies at z ~ 0.2. In this paper we have presented sev¬ 
eral good examples of its scientihc capabilities and potential. 
These results have been obtained using only a small fraction 
of the complete OMEGA dataset and therefore provide only 
a partial view of the star-format ion and AGN properties of 
the galaxies in the Abell 901/902 structure, biased towards 
the highest density regions. A more complete picture, cov¬ 
ering a much broader range of environments, will emerge 
once the full sample is analysed in future papers. Moreover, 
here we have only considered the integrated emission-line 
properties of the galaxies. OMEGA also provides spatially- 
resolved emission-line maps which will be used to study how 
the spatial distribution of the star formation is affected by 
internal and environmental factors. In future papers we will 
also explore the information on the chemical composition of 
the ionised gas encoded in the Ha/[Nil] emission-line ratio. 
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